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LINEAR EQUATIONS WITH TWO PARAMETERS* 



ANNA J. PELL 

The problem of the existence of solutions of two linear self -adjoint differential 
equations of the second order with two parameters has been treated by Hilbert, f 
Richardson, t and Yoshikawa.§ There are corresponding problems in the 
theory of linear equations in infinitely many unknowns, and in the theory of 
linear integral equations. Part I of this paper deals with the first problem, 
the existence of characteristic numbers X and ii for the system of equations 



Ui = X '^kij Uj + li^lijUj (t = l, 2,. . 0, 

(1) 

00 00 

1 = 1 1=1 

where the matrices are subject to certain conditions (§1), and the expansion 
of the determinant matrix in terms of the solutions. Part II deals with the 
second problem, the existence of characteristic numbers X and n for the equations 

u{x) = \ I K{x, y) u{y)dy + ix j L{x, y) u {y) dy, 

v{s) =\ JM {s, t) vit)dt- fjL J N{s, t) V (/) dt, 

where the kernels are subject to certain conditions (§3), and the expansion of 
arbitrary functions of two variables in terms of the characteristic functions. 

* Presented to the Society, October 25, 1919. 

t GrundzUge einer allgemeinen Theorie der linearen Integralgleichungen, pp. 262-267. 
X Theorems of oscillation for two linear differential equations of the second order with two 
parameters, these Transactions, vol. 13 (1912), pp. 22-34. 

§ Ein zweiparamelriges Oscillalionstheorem, Gottinger Nachrichten, 1910, 

pp. 586-594. 
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I. Linear equations in infinitely many unknowns 

1. Existence of solutions. Consider the system of linear equations (1) 
where the matrices K, L, M, N are real and symmetric, the sum of the squares 
of the elements in each matrix is convergent, the two matrices L and N are 
positive definite and feyM« + li/mki^O. We shall show that under these conditions 
there exist values of X and p. for which the equations (1) have solutions of finite 
norm and not identically zero, that is, neither [ui] nor {vk} is identically zero. 
The method consists of the elimination of /x from (1), the transformation of the 
resulting system of equations in the one parameter X into a system of the form 
(9) for which the existence of characteristic numbers X, has already been estab- 
lished, the substitution of these X, in one system of equations of (1), and the 
proof of the existence of corresponding characteristic numbers m^, and finally 
the proof that Xj and ju,a are also characteristic numbers for the other equations 
of (1). 

The elimination of m in (1) gives the system 



(2) 






Since L and N are positive definite their characteristic numbers are positive; 
call them cq and fi\. The matrices formed from the corresponding solutions 
l*j and nt, are orthogonal,* and hence the system (2) is equivalent to 

(3) uK (i + i)= xE(a=o| + -:< H\) vi 

where 

i t 

(4) K*=L*KL*' , M*=N*MN*', 

&ij = 1 if i = /, and = 0, if i^j. 
With the substitution 

/-■v 1 * Obi I * Otj 

(5) aj,=fey^+w.;,,-y, 

'^f 2 ^ ^' 



* For notation and terminology see Hilbert, loc. cit., Kap. XI, and Hellinger, E., Netie 
Begriindung der Theorie quadralischer Formen von unendlichvielen Verdnderlichen, Jour- 
nal fiir Mathematik, vol. 136, pp. 210-262. 
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the system (3) becomes 

(6) 
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2 * X '^"^ • 



The matrix A is symmetric, for the interchange of i with ; and of k with /, or the 
interchange of ^ with rj, leaves A unchanged. The matrix 

Cf a, + 0jfe a,- + 0k 

is completely continuous. Let X( = a;,,(, be of norm ^ 1 . Since 

<1. 



«. 



a? + -8? 



yj^ikkij 2 , o2^-'* = -1/7 Ai 2-(*'* 



and 



/ , ^iW^a 2 r o2 *^-'* 2-( ^■*'^0' 2 . o2 *i* 

"=1 ai+|3;i i.i,u=\ oci + Pk I 






By hypothesis S<y fe,;;' converges, and lim m-*8 S"j = „ k,j = 0. Hence the matrix 
k ij5ki oij/iat + 0k) is completely continuous, and similarly for mkfiijfik/ i.o^i + h)- 
This shows that a{,/C{ is completely continuous. From the inequality 



Hr, 



cf^ 



\Ct c./ 



it follows that a^^/cf^ is completely continuous. The matrix a^/cf^ is sym- 
metric, and not identically zero since KN + LM ^0, hence there exist* values 
of X for which the system 



(7) 



1 ''fv 



has solutions of finite norm. From this it follows that the system 
(8) z« = XE ^^' 



• Hilbert, loc. cit., p. 148. 
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has solutions [z^ = yf^} of finite norm for the same value of X. Since the ma- 
trix a(^/c( is completely continuous the adjoint* system 

(9) ^:=xz^< 

also has solutions {x*^} of finite norm. For each value of X there is only a 
finite number of linearly independent solutions of (8) and an equal number of 
linearly independent solutions of (9), and the relation between the solutions 
of (8) and (9) is 

(10) 0£ = 4x1 

A solution {x(} of (9) leads to a solution {xik } , such that S <_ ^ « j^ is con- 
vergent, of the system 

(11) Y^hjXjk + ^n kiXii = X ]^ {kijXjiHki + lijXjintki). 

We shall show later that Xn, either has the form UiV^, or is a sum of such terms. 
Return to the system (1), and let X be a value for which the system (11) 
has a solution Xn,. Assume first that M( is the only solution of 

(12) «■ = X ^ kijUj- 

J 

with the understanding that {u*} may be = 0, and that S,m*'' = 1 if m, is not 
= 0. There exists a limited matrixf P such that 



where 



P = \K + \PK - U*U*', 
w*fc=M* if k = l, Uik = if ^7^1. 



If the system 

(13) U{x)=\KU{x) + nLU{x) 

has a solution U{x) besides U*ix) for n = 0, then U{x) satisfies 

Uix) = ix{LUix) + PLU{x) ] + cU*ix), 

(14) 

{U*. LU)=0. 



* Hilbert, loc. cit., p. 165. 

t Hurwitz, W. A., On the pseudo-resolvent to the kernel of an integral equation, these Trans- 
actions, vol. 13 (1912), pp. 405-418. 
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Also if (14) has a solution Uix), this solution satisfies (13), and the two equations 

(13) and (14) have the same solutions except for U*{x). The two equations of 

(14) may be expressed by one equation, for from the second part it follows that 

c{U*, LU*)= - M[(f/*, LLU) + (U*, LPLU)], 
and the substitution of this value for c gives 

U*U*'LLU{x) + U*U*'LPLU{x)'l 



(15) U(x) = M llUix) + PLUix) 



([/*, LU*) J 

= l^QUix). 

The equations (14) and (15) have the same solutions except for U*{x). To 
establish the existence of characteristic numbers of (15), consider the system 



(16) C/(^) = m{ 



. LU*U* 'LU{x)+PLU*U*'LU{x) ' 
(J U (x) 



(U*, LU*) J 

= uRU{x). 

All the solutions of (15) satisfy (16), and any solution U(x) of (16) which is such 
that 

(17) {U*,LU) = 

satisfies (15). This condition is satisfied by all the solutions of (16) except possi- 
bly those corresponding to the value no of n given by 

^^[iU*, LLU*) + {U*, LPLU*)] + {U*,LU*) = 0. 

A solution of (16) corresponding to /no is U*(x), the solution of (13) for pi = 0. 
By adding multiples of U*(x) any other solution of (16) corresponding to juo 
can be made to satisfy (17), and hence we may assume that the systems (13) 
and (16) have the same solutions. The matrix R in (16) is such that the product 
of R by the symmetric positive definite matrix L is symmetric, and therefore* 
characteristic numbers m^ exist and are real, {LU^(x)} is the adjoint system 
of the system of solutions f/„(x) and the matrix R may be expressed in terms 
of the solutions in the following form : 

(18) R'LFix) = ^ LUjm^FJJ^)^ 



* A. J. Pell, Linear equations with unsymmetric systems oj coefficients, these Transac- 
tions, vol. 20 (1919), pp. 23-39. 
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where F{x) is any limited linear form. From this expansion it follows not only 
that (16) and (13) have solutions besides U*{x), but also that the system of all 
the solutions U„{x) of the equations 

(19) U^ix) =\KU„(,x} + p^J^U^ix) 

is such that the system {LU^ix) } is complete. Suppose that F{x) is such that 

(F,LUJ = 0; 

then from (16) and the expansion for R the linear form 

Gix)=LFix) + PLF{x) 
is such that 

Lf^n^K^^*)^^. 

{U*,LU*) 

The only G{x) which satisfies this isG'(a;) = kU*{x). From the definition of 
the matrix P it follows that if 

LF{x) +PLF{x)=kU*{x), 

then LFix) = cU*{x), and therefore Fix) = 0. Hence the system {LU^{x)} 
is complete. 

If the system (12) had more than one linearly independent solution, similar 
considerations would show that again the system {LU„{x) ] is a complete system. 

From (11) we obtain solutions V„ix) of the second system of equations of 
(1) by multiplying by C/„, 

(20) V„{x) =mVAx) - HoNV^ix), 

for the characteristic numbers X and jUa- The solutions V^ix) are given by 

V„{x}=X'LU„ixX 

and since {LU„{x)} is complete the Va{x) are not all zero, and the equations 
(19) and (20) have solutions not identically zero. We have established the 
following theorem. 

Theorem 1. If the matrices K, L, M, and N are symmetric matrices such that 
the sum of the squares of the elements in each is convergent, if L and N are positive 
definite, and if kijHki + hjWki ^ 0, there exist real values of X and n for which the 
system (1) has solutions m,-, v^ of finite norm and not identically zero. 
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2. Expansion of the determinant matrix. Since every solution m„,, v^k of 
finite norm of (19) and (20) satisfies (9), only a finite number of v^k can be 
different from zero. It has been shown that {LU^ix)} is a complete system, 
and therefore 

H 
a = l 

for otherwise the difference between the left and right hand sides would be 
orthogonal to LUa{x). We may therefore assume that the sohitionsof (11) 
are in the form of products Mai^ofc- 

From the expansion of the symmetric completely continuous matrix* (a^^/CfC,), 



it follows that 





a 




Of, = 


^2 


* 2 * 
Xa(C^X„^ 

K 



and this gives the following expansion for the determinant matrix: 

(21) k,jn,i + kptki = D ""^^f^, 

where 

i ' 

The following theorem has been established. 

Theorem 2. Under the conditions of Theorem 1, the determinant matrix KN 
+ LM may be expressed in terms oj the solutions of (1) in the form (21). 

II. Two LINEAR INTEGRAL EQUATIONS WITH TWO PARAMETERS 

3. Existence of solutions of two integral equations. In the two linear integral 
equations 



(22) 



M(ic) = X / K{x,y)uiy)dy + n I L{x,y)u{y)dy, 
v{s) = X / M{s,t)v{t)dt - n I Nis,i)vit)dt, 



* Hilbert, loc. cit., p. 148. 
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let the kernels K, L be real symmetric functions continuous in the real variables 
X and y for a ^ x ■^ b and a ^ y '^ b, M, Nreal, symmetric and continuous in 
the real variables 5 and t ior c ^ s ^ d, c ^ t ^ d, KN + LM ^ 0, and L and 
N positive definite in the sense that there exist no functions /(ar), g{s) not iden- 
tically zero which together with their squares are integrable in the sense of 
Lebesgue on (a, b) and (c, d) respectively, and such that J^ S f(x)L{x,y)f{y) 
dxdy ^ 0, f f gis)Nis,t)gU)dsdt ^ 0. Let { (p,(^) } be a closed orthogonal 
system of continuous functions for the interval (a, 6), and {i/'i(^)} for the in- 
terval (c, d). Multiply the equations (22) by <pi{x) and 4'k(.^)> respectively, inte- 
grate, and then expand the right hand sides; 

JikV = X Z J J 4'kMh J4'iv - fiY^ J j ^kNh J > 



(23) 



^tv. 



If the equations (22) have continuous solutions u{x), v{s), then the system (23) 
has solutions of finite norm. The matrices in (23) satisfy all the conditions 
imposed on the matrices in Theorem 1, and therefore there exist values of X and 
IX, necessarily real, for which the system (23) has solutions x„ y^ of finite norm. 
In the usual way it can be shown that the functions 



k J k J 



mk 



are continuous in (a, 6) and (c, d), respectively, and satisfy the equations 
(22). 

Theorem 3. // the kernels K, L, M, N are continuous real symmetric functions, 
L and N positive definite, and KN + LM^ 0, there exist values of \ and ix neces- 
sarily real and for which the equations (22) have continuous solutions u{x), v{s), 
not identically zero. 

4. Properties of the solutions. Let m,(«), Viis) be solutions of (22) corre- 
sponding to the parameter values X,-, m,-, and %(«), v^is) to X^, ma- Since the 
kernels are symmetric it follows from (22) that 
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/ UiUk = \ I I UiKuk + Hi I I UiLuk 

= ^k I I UiKuk -\- Hk I I UiLuk, 
(24) 

/ ViVk = \i I I ViMvk - Hi I I ViNvk 

= ^k I I VfMvk - Hk I I ViNvk- 
By subtraction we obtain 

{\ -h) I I UiKuk + (m. - M*) / / UiLuk = 0, 

(X,- - WjJviMvk - (Mi - Hk)JJviNvk = 0. 



(25) 



If the two parameter sets are not identical the determinant of the coefficients 
must equal zero: 



(26) 



/ / UjKuk I I UiLu 

J j ViMvk ~ f f ^•■^"* 



= 0. 



This relation together with (24) gives the result that if (m„ Vi), (%, Vk) belong 
to different parameter sets, the matrix 



(27) 



/ UiUk I I UiKuk I I Ui 

I 1>iVk I I ViMVk ~ I I '"i' 



is of rank less than 2. If the solutions belong to the same values of the 
parameters, linear combinations may be formed so that this condition is 
satisfied. 

Since the kernels L and N are positive definite the solutions Mj, b,- may be 
multiplied by constants such that 

/ mW / ViNvi + j v'^i j j UiLui = 1. 
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In the following we shall assume that the solutions have been so normal 
ized. 

Consider any set of functions «,(«), Vi{s) which are continuous on (a, b), 
(c, d), respectively, and have the property that 

(28) / UiUk I I ViNvk + J VfVk I I UiLuk = 5ik- 

If the series 

(29) J{x,s) = Y^Aflii{x)vi{s) 

i 

is uniformly convergent the coefficients Ai may be determined in terms of 
f{x, s). Multiply the equality (29) by 

Tk{s) I L{x, y)uk{y)dy + Uk{x) I N{s, t)vk{t)dt 

and integrate; then on account of (28) . 

Ak= I I I Uk{x)L{x,y)f{y,s)Vk{s)dxdyds 

+ I I I Ukix)f{x,s)N{s.t)uk{t)dxdsdt. 
Similarly if the series 

fix.s) = I^5.h(5)/'«,e>')L(>',;c)dy + M.(^)yiV(5.<)»<(0c?«] 
is uniformly convergent, 

Bi= j I Ui{x)f(x,s)vi{s)dxds. 

Let /(a;, s) be any function continuous in x on (a, b) and in s on (c, d), and let 

/.= f f fui{x)L{x,y)f(y,s)vi{s)dxdyds + J J J Ui{x)f{x,s)N{s,t)vi{t)dxdsdt. 

Call/, the Fourier coefficients of /(^, s) with respect to %, w,. If (29) is uniformly 
convergent, 

Jjigi=ffff{x,s)L{x,y)g{y,s)dxdyds+JjJg{x.s)N(s,mx,t)dxdsdt, 
where g{x, s) is a continuous function. 
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The Fourier coefficients of any continuous function f(x, s) with respect to a 
set of functions m,-, Vi with the property (28) are of finite norm. For 

(30) / / I F{x,s)L{x,y)F{y,s)dxdyds + J I j F{x,s)N{s.t)F{x,t)dxdsdt^Q, 
and if 

n 

F{x,s)=J{x,s) - ^fiUi{x)vi{s), 
the inequality (30) reduces to 

J I j f{x,s)L{x,y)f{y,s)dxdyds + f J f Ax,s)N{s,t)f{x,t)dxdsdt 

Therefore 

(31) X'fi^J J Jf{x,s)L{x,y)f{y,s)dxdyds ^ J J Jj{x,s)N{s,t)j{x,t)dxdsdt, 

which shows that the sequence {/^ } is of finite norm. 

5. Expansion of arbitrary functions of two variables. A consequence of 
(23) and (21) is 

(32) f J f f{K{x,y)f{x,s)N{s,t)g(y,t) + L(x,y)f{x,s)M{s,t)g(y.t) jdxdydsdt 



a "a 



where /(a;, s) and g(x, s) are any two continuous functions, and therefore if the 
series on the right is uniformly convergent, 

Kix,y)N(s.t) +L(^,y)M(5,0=i: ^''^^'f^°^'^\ 
where 

W„ix,s)=v„is) j L{x,y)u„{y)dy + u^{x) j N{s,t)vMdt. 

To obtain the expansion of arbitrary functions, we first show that 
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is of finite norm; this would follow immediately if there were a continuous func- 
tion F{x, y, s, i) such that 

III I Hx,y,s,t)<Pi{x)<pj(y)\f>kis)Mt)dxdydsdt = bikji = aikji-^^-^, 
aiiu = f f K{x,y).Pi{x).p^{.y)dxdy ^' + J J M{s,t)Us)h{t)dsdt% 



where 



and <Pi{x) are the characteristic functions of L{x,y) corresponding to the charac- 
teristic numbers a?, and ^ki.^) the characteristic functions of N{s, t) corre- 
sponding to the characteristic numbers /3^, for then 

WAx,s)=Kf f P{x,y,s,t)W„{y,t)dydt, 

and (31) would give the convergence of 

^ Wl{x,s) 

Denote by F\J{x, s)] the transformation defined by 

/ I F[f(x,s)],pi(x)4'k(s)dxds = Ylbikji I j f{x,s)vj{x)^i{s)dxds. 
If the function F{x, y, s, t) above exists, then 

F[f{x, s) ] =JJP('^' y-^' 0/0'. t)dydt. 

The transformed function of g{s)J'L{x, y) f{y)dy exists and is a continuous 
function, for it may be expressed by 

// Kix, y)pi(y)dy I f<pi 2 /• 

L{x,y)f(y)dy] = X) ^ 5 ^^Ms) / g4^k 

<Pi(x) I ffi/ fii r r \ 

+ 1: 4— (-yf^ / g^J M{s.t)Mt)dt\ 

Xk on Xoi+PkJ J f 

a series with continuous terms and such that the series formed from the absolute 
values of the terms is uniformly convergent.* Similarly the transf ormed f unc- 

* Mercer, J., Functions of positive and negative type. Transactions of the 
Royal Society, A, vol. 209 (1909), pp. 415-446. 
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tion oif{x)J^N{s, t) g{t) dt exists and is continuous. From the previous work 
it follows that the transformed function of W^ix, s) is given by 

(33) F[W^{x,s)]= ^"^^''^ 



K 

From the series which represents F[g(j)y"L/] the following inequality is obtained : 



and hence 



(34) 2\F{gjLf]\^JjLjJg' + ^Kx,s), 
where 

iKx,s) = 2Y, ( fKix,y).p,(y)dyy-p^+^ ( f MisMkiOdt)' 
i.k \»^ / ftj + Pk at \J / 

and is a continuous function in x and s. Similarly it may be shown that 

(35) 2\F\fjNg] 1 ^ fs^&fP + *'(^. ^). 

where •^'^{x, s) is a continuous function in x and s. From the two inequalities 
(34) and (35) it follows that 

*X^, s) - 2F[J L(x, y)f(y. s)dy] +ffff(^' ^)^(^. 3')/(5'. s)dxdyds 
(36) 

+ 'ifKx.s) -2F[jN{s,t)fix,t)dt] + ffff(y' ^)^(^- ')/(?'• t)dydsdt ^ 

where „ 

Ax, s)=Y^F/Uj{x)vj{s) 
J-i 
and 

or by virtue of (33) 

Fj=FIvj(s)Jl(x, y)uj(y)dy + uj{x)J'n{s, t)vj(t)dt]. 
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The inequality (36) reduces to 

J-i j=i 

and therefore 



This last result and the inequality (31) give the uniform convergence of the 
series formed from the absolute values of the terms of the series 

where f{x, s) is any continuous function, and from (32) it follows that 

h{x, s) =ff{K{x, y)J{y, t)N(s, t) + L(x, y)f(y, t)Mis. t))dydt 
^Y^W^^^Y.WAx,s)fMj,.. 

Theorem 4. // the kernels K, L, M, N are continuous real symmetric functions, 
L and N positive definite in the sense defined in §3, and K N -\- L M ^Q, any 
function h{x, s) which can be expressed in the form 

h{x, s)=ff{K{x, y)f(y, t)Nis, t) + L{x, y)f(y, t)Mis, t))dydt, 

where f{x, s) is any continuous function, may be expanded into the uniformly and 
absolutely convergent series 



h(x, 5) = 2^ W„ix, s) I huj)„ 



Bryn Mawr Cohege, 
Bryn Mawr, Pa. 



